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Abstract 

The magnetism of superconductors has not been fully explored despite many years of research 
and its great significance to both theory and applications. Here we present a unified perspective on 
how quantum electromagnetic vortices and electron-spin textures join to protect superconductivity 
against the destructive power of large magnetic fields. Using rapidly rotating ultra-cold gasses of 
fermionic atoms as a model system amenable to detailed experimental exploration, we discover a 
hierarchy of Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phases in which a metal or a band-insulator 
of unpaired particles coexists with a spatially modulated superfluid. Quantum fluctuations can 
transform these phases into strongly correlated vortex liquid metals or insulators. We elucidate a 
general criterion for the stability of FFLO phases and discuss implications for the quantum Hall 
and vortex liquid phenomena in the pseudogap state of high temperature superconductors. 
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Some fundamental aspects of magnetism in superconductors remain unclear despite many 
years of research. Several mechanisms for the survival of superconductivity in large magnetic 
fields have been studied: vortex lattice formation [1], re-entrant superconductivity {2], J, 4], 
and the relaxation of Zeeman effect through spontaneous development of spatial modulations 
FFLO states) jsl, Q]. Only the first mechanism has been unambiguously observed in nature 
3] . Fortunately, a new frontier has been opened in atomic physics to observe and study these 
phenomena. Here we analyze in detail the very rich phase diagram of rotating fermionic 
atoms with attractive interactions near unitarity in two dimensions. ^ 

Electronic systems are very similar to trapped ultra-cold fermionic alkali atoms [8]. The 
orbital effect of magnetic field can be created in atomic systems by rotation. Coriolis forces in 
the rotating frame implement a uniform flux density which can produce quantized vortices 
and quantum Hall effect. The paramagnetic or Zeeman effect is created by preparing a 
system with unequal number of atoms in the two lowest-energy hyperfine states, which 
essentially have the same dynamics as electron spin and will be referred to as "spin" from 
now on. The Zeeman and orbital effects can be independently controlled. The strength 
of attractive interactions between atoms, responsible for superfluidity, is easily tuned using 
Feshbach resonances. Therefore, cold atoms are ideal systems for exploring a wide range of 
phenomena related to superconductivity. 

The competition between superfluidity and Zeeman effect is normally resolved by form- 
ing either an unpolarized superfluid of Cooper pairs, or a spin-polarized normal state 



The elusive FFLO state is a compromise in which the superfluid and spin- 



polarization coexist but alternate in space in a periodic pattern. FFLO states have been 



pdependently studied .n conde^edmatterandhidi energy physics |13,. Three-d„ne„si„na, 



|l4l . Il5| and two-dimensional Q, [l^, Q, Q, l2Q|, l2j| FFLO states have been considered in 



superconductors and heavy-fermion metals. The latter studies were motivated by the lay- 
ered structure of materials, which allows controlling independently the orbital and Zeeman 
effects by tilting the magnetic field with respect to the layers. In order to directly address the 
role of these effects in stabilizing two-dimensional FFLO states, we depart from the earlier 
studies and focus on the lowest temperatures where Landau level quantization is important. 
This will also establish a connection between re-entrant superfluidity, FFLO phenomena 
and strongly correlated insulators of general interest in unconventional superconductivity, 
using a previously unexplored rotating cold-atom perspective (for a review of past research. 
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mostly on rotating bosons or three-dimensional fermionic systems, see [22] and references 
therein) . 

Here we discover a sequence of FFLO phases at zero temperature in which a metal or a 
band insulator of unpaired atoms coexists with a superfluid hosting a vortex lattice. Sur- 
rounding these phases are strongly correlated vortex liquid states and atomic quantum Hall 
insulators. Our results clarify the most general stabilizing mechanism for FFLO states and 
provide a unified perspective on the fate of superconductors in high magnetic fields. We 
also elucidate the potential of trapped rotating ultra-cold atoms to be the first system that 
clearly realizes the FFLO states. This and other studies in the atomic physics setting are 
likely to shed light on the possible role of FFLO states in superconductors, heavy- fermion 
metals, quantum chromodynamics and astrophysics. Finally, our results point to experi- 
mental circumstances for the cold-atom realization and exploration of strongly correlated 
states which shape the unconventional properties of high-temperature superconductors. 

The quantitative applicability of the following universal theory to ultra-cold fermionic 
atoms near the Feshbach resonance stems from the unitarity limit. The scattering length 
of two-body interactions diverges at the resonance, so that the microscopic details of in- 
teractions become irrelevant and all properties of the many-body system become functions 
of only dimensionless ratios of externally controlled parameters. Such scaling laws are pro- 
tected inside all superfiuid phases, but may break down due to rotation in two dimensions 
when strongly correlated vortex liquid states are formed [23] . 

We consider a Bogoliubov-de Gennes (BdG) model in grand-canonical ensemble which 
arises as a saddle-point approximation in the quantum field theory of neutral fermionic par- 
ticles near unitarity [23! . This approximation becomes quantitatively accurate at sufficiently 
large densities or sufficiently far from the Feshbach resonance in the BCS (Bardeen, Cooper, 
Schrieffer) limit, but remains extremely useful even in the numerically accessible regimes for 
elucidating various qualitative features which may be observed in experiments. Rotation at 
angular velocity uj is described by a static gauge field A in the rotating frame which satisfies 
V X A = 2mujz^ where m is the mass of atoms. We assume that the centrifugal forces are 
almost exactly cancelled by the potential which traps atoms. The Zeeman effect occurs in 
finite Zeeman field h which together with the global chemical potential /x sets the chemical 
potentials for the two hyperfine species, /a^ = /a + h and jiy = ji — h. 

In the absence of interactions all atoms would be localized in their Landau orbitals. Such 
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localization facilitates Cooper pairing and allows superfluidity to occur even in the deep 
quantum limit at high rotation rates [2|, |3|, IJ]. The saddle-point approximation regards the 
superfluid as a macroscopic condensate of Cooper pairs into a single state given by the order 
parameter function: 

The amplitudes 0^^/ describe the occupation of bosonic Landau levels n G {0, 1, 2 ... } at 
wavevectors = integer x Sq in x-direction. The Landau level wavefunctions of dimension- 
less coordinates ^ = {^x^Cy) 

ry) = ^i=e^^«-e-^(«^+'')'//„(4 + v) (2) 

in Landau gauge, where are Hermite polynomials. This realizes a vortex lattice with 
periods Ax = 27Th/6q and Ay = ny6q{4:muj)~^ containing riy vortices per unit-cell. The 
hexagonal Abrikosov lattice has Uy = 2^ 5q = {4:7rV3mh(jjy^'^ ^ and (j)n^i = {—l)^i(j)n,o- 

The BdC Schrodinger equation determines the quantum-mechanical spectrum of 
fermionic quasiparticle excitations in the condensate: 
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Here, and are the particle and hole wavefunctions of opposite spins in the eigenstate 
with quantum numbers n. We solve this equation in the basis of Landau orbitals using 
exact diagonalization. The energy spectrum determines the free energy J-{(j)n^i) in the 
saddle-point approximation |i24l |: 



^{Ki) = I'^bTIu (1 + e-^(^"(^-^)+^)) + fc^rin (l + e-^(^"(^-^)-^)) (4) 



+ 



no,o(o,o) 



^=h=T=0 



dV|$(r)|2 . 



Here, (3 = {kBT)~^^ ks is Boltzmann constant, and u = —a^/a measures the strength of 
attractive interactions between particles through the two-body scattering length a and the 
conflnement scale in z-direction which renders the system two-dimensional. Therefore, u is 
proportional to the detuning from the Feshbach resonance in experiments. The ultra-violet 
behavior is regulated by rotation, and the subtracted pairing susceptibility Iln^^n^{iuj^Px) 



assures that u = corresponds to unitarity [23] . The global minimum of J-'{(/)n,i) determines 
in the thermodynamic equilibrium. Note that this approach is not limited to the vicinity 
of second order transitions like the earlier studies. 

The BdG spectrum E^i in the superfluid at /i = has a particle-hole symmetric band- 
structure reminiscent of broadened Landau levels at high energies. The low energy bands 
can be associated to quantum tunnelling of quasiparticles between vortex cores. The Zeeman 
field h acts as a chemical potential for the BdG quasiparticles which breaks the particle-hole 
symmetry. At zero temperature all states with < h are occupied, and with > h 
are empty. There is a "metal-insulator" quantum phase transition each time h crosses a 
band-edge [25]. As long as superfluidity survives, the extra quasiparticles at energies < 
En < h carry spin polarization which is mostly concentrated near vortex cores. Translational 
and rotational symmetries are spontaneously broken by both the superfluid density (vortex 
lattice) and the coexisting spin-polarization, therefore these states are FFLO states. 

In the traditional FFLO picture the order parameter spontaneously develops spatial mod- 
ulations in response to Zeeman effect alone. If the lowest positive energy in the resulting 
quasiparticle band-structure is larger than the pairing gap, then the FFLO state cannot be 
stable because it would breach the Pauli-Clogston limit, making it more advantageous to 
align all spins with the Zeeman fleld than to form singlet Cooper pairs. This problem is 
avoided in the presence of a vortex lattice. The energy gap for the vortex core states is 
always lower than the pairing gap (see FigH]). 

We numerically minimized the free energy density (|4j) using a fermionic basis of 72 
quasi-periodic states in each of the 50 lowest Landau levels. The global minimum was 
calculated among order parameters within the lowest flve bosonic Landau levels and = 2, 
5q = {4:7^^/3mh(JJy/'^^ allowing a variety of lattice structures. All results are summarized by 
the zero-temperature phase diagram in FigM 

The superfluid is divided into a number of phases with either gradually varying polar- 
ization, or polarization saturated at an integer number of spins per vortex. The unpaired 
polarized atoms form a Fermi surface in the former case, and a band-insulator in the lat- 
ter case. At least some of the "metal-insulator" transitions are flrst order in terms of the 
superfluid order parameter, which allows them to persist at flnite temperatures despite the 
loss of the sharp distinction between "metals" and "insulators" of the vortex core fermions. 
The phase diagram contains a sequence of quantum tri-critical points where the transitions 
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between atomic quantum Hall insulators meet the first and second order superfluid-insulator 
transitions. The superfluid boundary retains its pair-breaking first order transitions below 
some finite temperatures. 

Interestingly, Zeeman effect can produce FFLO superfiuids condensed in higher bosonic 
Landau levels [21]. Typically, one Landau level is much more populated than the other 
ones, as shown in Figl2] (a) and (c). Higher Landau levels can become favored at large 
h because they introduce additional vortex-antivortex pairs and therefore have larger re- 
gions of suppressed superfiuidity where local spin polarization can be significant (see Figl3]). 
Structural transitions of the vortex lattice within the same Landau level are also found in 
the saddle point approximation, but quantum fiuctuations of the order parameter generate 
familiar logarithmic "Coulomb" forces between vortices which generally act to stabilize the 
hexagonal Abrikosov lattice if it occurs in the lowest Landau level. 

Finally, re-entrant superfiuidity is refiected in the robust paired states which occur when- 
ever the chemical potentials of the two atom species approach Landau levels. This is best 
demonstrated in Figure [4] which shows pairing free energy. As a result, a sequence of paired- 
state islands extends to arbitrary large h and survives at /i = regardless of how weak 
the attractive interactions between atoms are. The experimental consequences are discussed 
below. 

The limitations of saddle point approximation (SPA) used to calculate the phase diagram 
in Figl2] are discussed in [23|]. An indicator for the quantitative applicability of SPA is 



n^^ ^ 1, where n is atom density and ^ ^ ^2/x/mA2 the BCS coherence length (A is 
the spatially averaged BCS gap). This quantity grows rapidly with density or detuning as 
shown in FiglHl The only qualitative shortcoming of SPA are the missed vortex liquid phases 
which are strongly correlated insulators of Cooper pairs with non-universal properties. These 
phases inevitably reside in the vicinity of all superfiuid-insulator transitions shown in Figl2] 
and grow at the expense of the vortex lattice regions, especially at the shown low densities 
when the number of atoms is not much larger than the number of vortices. The SPA still 
qualitatively captures the boundary between the paired correlated insulators and unpaired 
atomic quantum Hall states, as well as the bosonic Landau level occupation, and indirectly 
predicts the existence of exotic metallic and insulating spin-polarized vortex liquids. The 
first order transitions between vortex liquids at different densities and low temperatures are 
also expected to descend from the T = "metal-insulator" transitions. We leave for future 
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research the more detailed studies of vortex hquids, which are potentiaUy responsible for the 
unconventional properties of strongly correlated states in high-temperature superconductors 



26]. 



The temperatures needed to explore the vortex-FFLO phase diagram in cold atom ex- 
periments should be below T* = k^^hu^ which is a few nK for typical u ^ 27t x 100 Hz 
and within reach by evaporative cooling. The main difficulty would be working with a low 
enough density of trapped atoms, ideally not much larger than no = muj/{T{h) ~ 10^ — 10'' 



cm ^, the density for a single populated Landau level [27|,l28|]. However, the main features of 



the phase diagram in Figl2] are expected to persist at much higher densities than shown, so 
that achieving the deep quantum limit may not be necessary. New experimental techniques 
may be able to overcome this difficulty in the future {29I. Figl6] illustrates some possi- 
ble experimental situations estimated from the grand-canonical results using local density 
approximation (LDA). 

In this paper we analyzed a generic class of stable FFLO states in two-dimensions at zero 
temperature, whose inffuence is observable at finite temperatures. These states result from 
the combination of orbital and Zeeman effects which can be independently controlled in cold 
atom experiments. Further research, especially cold atom experiments near unitarity, could 
provide valuable answers to long standing questions about the stability of superconducting 
phases in large magnetic fields and the nature of strongly correlated normal states. 

The mechanism for FFLO states described in this letter is perhaps the most generic 
and robust one. It requires a gap in the Bogoliubov quasiparticle spectrum smaller than the 
superfiuid pairing scale. A vortex lattice is a natural agent to realize this situation due to the 
existence of sub-gap core states. One question is whether spontaneously developed spatial 
modulations in time-reversal symmetric conditions can produce such a quasiparticle band- 
structure that the band gap at the chemical potential is smaller than the superconducting 
gap. The answer to this question could be important for finding a stable traditional FFLO 
state. 

I am very grateful to R.G.Hulet, Z.Hadzibabic and L.Baksmaty for useful comments and 
discussions. Numerical calculations were performed on Rice University computing clusters. 
This research was supported by W.M.Keck Program in Quantum Materials. 
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FIG. 1: Vortex core bands in relation to the superfluid gap near a vortex. The quasiparticle 
spectrum reduces to nearly-degenerate "Landau levels" at high energies, but resembles a band- 
structure shaped by the vortex lattice periodic potential at low energies, especially below the gap 
scale Aq. 
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FIG. 2: Quantum phase diagrams (T = 0) of fast-rotating fermionic atoms near unitarity, showing 
spatiahy averaged density of Cooper pairs ps and magnetization m — — ny. (a) ps and (b) m at 
unitarity; (c) ps and (d) m in the BCS hmit dX v — 0.15. The density plots of ps are color-coded by 
the dominant Landau level (LL) n in the condensate, also indicated with numbers: red n = 0, green 
n — 1. The contour plots of m emphasize superfluid-insulator transitions with yellow lines (straight 
vertical segments are 2^^ order, curved segments are 1^* order), metal-insulator transitions within 
the superfluid phase with red lines (dotted in the ps plots), and transitions between quantum 
Hall insulators (I) with thick-black lines. The superfluid phases (SF) host an Abrikosov vortex 
lattice, and may contain a coexisting metal (SF-M) or band-insulator (SF-I) of polarized unpaired 
fermions. The numbers in the phase labels indicate the highest populated quasiparticle band. Note 
that "spin-density" is saturated at an integer number of particles per flux quantum in all SF, SF-I 
and I phases, while it is not quantized in SF-M phases. 
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(a) (b) (c) (d) (e) 




FIG. 3: Density plots of hexagonal vortex lattices for condensates in (a) n = 0, (b) n = 1, (c) 
n = 2, (d) n = 3, and (e) n = 4 Landau level. Brightness is normalized separately in each plot. 
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(a) (b) 

FIG. 6: Radial (r) profiles of Cooper pair density Ps(r), total particle density n(r) and polarization 
p{r) as a function of total trap polarization P = {N^ — Ni)/N. The total particle number N = 
+ in the trap is = 3.5 • 10^; dimensionless detuning is (a) = 0, (b) = 0.15. The 
total polarization P is indicated on the top. These plots illustrate sharp density features at very 
low temperatures. Starting from P = the initial trend is that the superfluid core shrinks as the 
added "majority-spin" atoms are pushed to the trap boundary. However, the ring structure quickly 
acquires complexity from the FigE moving from r = toward the trap edge follows a trajectory of 
reducing fi while keeping h constant in FiglJl Characteristic superfluid rings separated by normal 
regions are a signature of re-entrant superfluidity, and polarization can be finite in a superfluid 
ring (FFLO states, or SF-I and SF-M in FigEj). The distinction between SF-I and SF-M phases 
can be made by wheather p{r) is a quantized constant or smoothly varying respectively. The SF 
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regions replaced by vortex liquids are expected to have strong pairing fluctuations [30|, 
currently experimentally unreachable trap harmonic frequency (jJ± = (1 + 5 • 10~^)(jJ was chosen to 



justify local density approximation, and the shown spatial region is r < 500 \/h /(2m(jj). All other 



parameters are comparable with already accomplished experimental setups 
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